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Abstract 



I show that the cluster variation method, long used as a powerful hierarchy of 
approximations for discrete (Ising-like) two-dimensional lattice models, yields 
exact results on the disorder varieties which appear when competitive interac- 
tions are put into these models. I consider, as an example, the plaquette ap- 
proximation of the cluster variation method for the square lattice Ising model 
with nearest-neighbor, next-nearest-neighbor and plaquette interactions, and, 
after rederiving known results, report simple closed-form expressions for the 
pair and plaquette correlation functions. 
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I. INTRODUCTION 



The cluster variation method (CVM) is a powerful hierarchy of approximations for lattice 
models of equilibrium statistical mechanics which has been invented by Kikuchi |IJ and more 
recently rewritten by An and Morita || . It is particularly well suited to analyse complex 
phase diagrams of discrete classical models [|J, but in some simple cases it is also known 
to give exact results. Since the approximations involved amount to neglecting correlations 
except for a finite range, exact results are obtained whenever correlations have a particularly 
simple structure, as in tree- like lattices f^-f| or one- dimensional strips |J. 

The purpose of the present paper is to study the behaviour of the CVM in another 
situation in which correlations are particularly simple, namely in the case of disorder varieties 
of two-dimensional Ising-like models with competitive interactions. Disorder varieties are 
known since the papers by Stephenson [0)0] and have subsequently been studied by many 
authors ]T^-|TU| . On a disorder variety (which is a suitable subspace in the whole parameter 
space of a model) the correlation functions factorize in a simple way, which leads to an 
effective dimensional reduction of the model, so one could expect that the CVM might be 
particularly accurate or even exact in such a case. This is indeed the case and I shall show, 
giving both general arguments and a detailed analysis of a particular model, that the CVM 
is exact on disorder varieties. 

The plan of the paper is as follows: in Sec. II I shall introduce disorder varieties and 
briefly recall some of the results which have been obtained in the past years; Sec. Ill will be 
devoted to the definition and explanation of the CVM; in Sec. IV the exactness of the CVM 
on disorder varieties will be shown and finally, conclusions will be drawn in Sec. V. 



II. DISORDER VARIETIES 

A disorder variety is a subspace of the parameter space of a model with competitive 
interactions, lying in the disordered phase, where the correlations have a particularly simple 
form and the model can then be integrated exactly. The first example of such a variety 
has been found by Stephenson |Tl| in the anisotropic antiferromagnetic Ising model on the 



triangular lattice. The hamiltonian of the model can be written in the form 

H = - Y,JijVi<Tj, (1) 

(ij) 

where <7j = ±1 is the spin variable at site i, the sum is over all nearest-neighbour (NN) pairs 
and Jjj depends only on the direction of the link between sites i and j. The values of Jjj 
along the three lattice directions will be denoted by Ji, J2 and J3. In the antiferromagnetic 
model we have J/ < for / = 1, 2, 3. Stephenson showed that when the condition 

tanhA^ + tanhA'xtanhA'a = 0, K x = Ji/k B T (2) 

(or one which is obtained from it by a cyclic permutation of the indices) holds, then the pair 
correlation along a lattice direction has a simple exponential form, as for the one-dimensional 
model. If Oi and aj are two spin variables separated by a distance fcona linear chain of the 
lattice in the Zth direction, their correlation ((Ji(Jj) is the fcth power of the NN correlation 
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along the same direction. In particular, assuming J\ < J 2 < J3 < 0, one has (04 cr.,) = 
[tanh(if 1 )] fe in direction 1, (<Ji<Jj) = [tanh(i^ 2 )] A ' in direction 2 and ((TiO-j) = [tanh(— i^ 3 )] fc 
in direction 3. Stephenson also showed that the disorder variety separates a portion of the 
disordered phase in which the pair correlation has an oscillating behavior from one in which 



it decreases monotonically. Similar results have been obtained by the same author |TT| for 
the union jack lattice and for certain one-dimensional lattices. 

Later, Enting |E| showed that the interaction round a face (IRF) model on the square 
lattice (and in particular the Ising model with NN, next-nearest-neighbour (NNN) and pla- 
quette interactions) has a disorder variety which can be mapped onto an exactly solvable 
crystal growth model. Peschel and Emery |]13[ rederived Stephenson's results for the corre- 
lations on the disorder variety of the triangular Ising model by means of a one-dimensional 
kinetic model and applied this technique also to the ANNNI model. Peschel and Rys |14j 
solved the eight vertex model on one of its disorder varieties. 

Baxter [JTSJ analysed the disorder varieties of the IRF model on the square lattice. He 
showed that the eigenvector of the (diagonal to diagonal) transfer matrix corresponding to 
the largest eigenvalue can be written in a simple form as the product of a sequence of two-site 
(NN) factors. 

Rujan |Ll| studied the relations between different techniques and considered several mod- 
els (vertex models, staggered IRF model, g-state Potts models, random bond models). 

Jaekel and Maillard |T/J found a local criterion which characterizes disorder varieties for 
any dimensionality and explains the effective dimensional reduction occurring in the model: 
the Boltzmann weight of an elementary cell of the lattice, summed over some (suitably cho- 
sen) spins (or whatever degrees of freedom), is independent of the remaining spins. Georges 



and coworkers 18] used this local criterion to calculate correlation functions on the disorder 



varieties of three-dimensional Ising models. 

To conclude this (certainly not exhaustive) brief survey of the existing literature, we 
mention that recently, Meyer and coworkers studied the disorder varieties of the eight 
vertex model in the framework of a random matrix theory approach to the transfer matrix. 



III. THE CLUSTER VARIATION METHOD 

The cluster variation method (CVM) is a hierarchy of approximation techniques for 
discrete classical lattice models, which has been invented by Kikuchi 0. In its modern 
formulation [plRl the CVM is based on the truncation of the cumulant expansion of the 
variational principle of equilibrium statistical mechanics, which says that the free energy T 
of a model defined on the lattice A is given by 

T = min F[ph\ = min Tr (p A H + p A lnp A ), (3) 

where H is the hamiltonian of the model, k^T = 1 for simplicity, and the minimization 
must be performed with respect to a density matrix obeying the normalization constraint 
Tr(p A ) 1. 

If the model under consideration has only short range interactions and the maximal 
clusters are sufficiently large the hamiltonian can be decomposed into a sum of cluster 
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contributions H a and the approximate variational free energy takes the form 

F[{p a , aeM}]=J2 [MpM - a a S a ] , (4) 

where a is a cluster of sites, p a = Tr A \ a pA is the cluster density matrix (TrA\a denotes 
a summation over all degrees of freedom except those belonging to the cluster a), S a = 
— Tr(p Q ha. p a ) is the cluster entropy and the coefficients a a can be easily obtained from the 
set of linear equations [0,|3| 

£ a a = l, V/3 6M. (5) 

f3CaeM 

The cluster density matrices must satisfy the following normalization and compatibility 
conditions 

Trp a = l,Va G M and p a = Trp^pp, Va C (5 6 M. (6) 

Notice that (f|) would still be exact if the density matrix pa of the whole lattice could 
be written exactly as a product of cluster density matrices in the form 

Pa = IW (7) 



IV. EXACTNESS OF THE CLUSTER VARIATION METHOD ON DISORDER 

VARIETIES 

There are two properties of the disorder varieties which suggest, at least for two- 
dimensional models, that the CVM might be exact on them. One is the one-dimensional-like 
character of the pair correlations. In fact, it is known that for a one-dimensional model with 
NN interactions, the pair approximation of the CVM (that is, the approximation in which the 
maximal clusters are the NN pairs), which is equivalent to the Bethe-Peierls approximation, 
is exact |],|20 . 



The other property, still valid for two-dimensional models, is related to a result by Baxter 
15fl . He showed that the eigenvector (corresponding to the largest eigenvalue) of the diagonal 
to diagonal transfer matrix is simply the product of a sequence of two-site (NN) factors. Since 
the density matrix of a diagonal cluster is the square of this eigenvector also the density 
matrix has a product structure. As we have seen in the previous section, when the density 
matrix has a suitable product structure the CVM becomes exact. Therefore one can hope to 
find a CVM approximation which is exact on the disorder variety of a given two-dimensional 
model. In the square lattice case a good candidate is the plaquette approximation |1,21 



which is equivalent to the Kramers- Wannier approximation |f22fl , which in turn has long been 



known to correspond to a variational approximation in which the largest eigenvalue of the 



transfer matrix is searched using a restricted space of factorized vectors |23 . 

In the case of the plaquette approximation for a model defined on the square lattice the 
condition (|7|), which implies the exactness of the approximation, becomes 

p Oplaq Pplaq Ilsite Psite 

Ilpair Ppair 
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where p\ denotes the density matrix of the whole lattice and the products are to be intended 
over all plaquettes, pairs and sites of the lattice. The above equation should however be 
taken with some care, since it is known that not all local thermodynamic states (i.e. density 



matrices) can be extended to the whole lattice p4j. Consider as an example a model of 
Ising spins Oi = ±1 in its disordered phase, which will be studied in detail below. One can 
easily check on small lattices that, using a generic plaquette density matrix and the pair 
and site matrices derived from it by partial traces, (H) leads to a pa which is not correctly 
normalized. In the case of open boundary conditions (with this choice the sites and the pairs 
lying at the boundary do not enter the products in (^)) the correct normalization is achieved 
only if d = c 2 , where c = ((XjCr^NN an d d = (c^cr,- )nnn are the NN and NNN correlations, 
respectively. 

When the condition d = c 2 holds, the procedure of extending local density matrices 
to larger clusters is well defined, in the sense that by partial traces one can reobtain the 
local density matrices which were used to build the larger ones. In addition, one can verify 
that the density matrix of any cluster admits a decomposition into a product of plaquette, 
pair and site density matrices, with exponents given by the CVM rules. For instance, with 
reference to Fig. [l|, in the case of the 3x3 square we have 

/ x _ Ppl a q(Tl,T2,T5,T4)ppi aq (T 2 ,T3,T6,T5)ppi aq (T4,T5,T 8 ,T7)ppi aq (T 5 ,T6,T9,T 8 )p s ite(T 5 ) 
Pq{Ti, . . .Tg) — 



Ppair(T2, T 5 )p pa i r (t 5 , T 8 )p pair (t 4 , T 5 )p pair (T 5 , T 6 ) 



(9) 



while for the zig-zag chain 



/ N Ppair(0"l,0"2)Ppair(0"2,0"3) • • •Ppair(0"L-l ) 0"i) , ln x 
Pchain (01,0-2, . . .a L ) = r r - r . (1U) 

Psite(C r 2)Psite(,0'3j • • • Psite {&L-1 ) 

As a consequence, also the pair correlation function has a very simple product form, that is 
(labeling the spin variables by the site coordinates) 

g(x, y) = (a(x Q , y )cr(x + x,y + y)) = c N+|y| . (11) 



The result ([U]) is equivalent to the result by Baxter |T3| that, on disorder varieties, 
the eigenvector of the diagonal to diagonal transfer matrix, corresponding to the largest 
eigenvalue, can be written as a product of NN pair terms. p c ham is just the square of this 
eigenvector, and the site factors which appear in the denominator can be easily associated, in 
a symmetric way, to the adjacent pairs. This shows (although this is not a rigorous proof) 
that when the plaquette approximation is exact for a model of Ising spins on the square 
lattice, then the model is at a point of the disorder variety in its parameter space. 

Let us finally study in detail the square lattice Ising model with NN, NNN and plaquette 
interactions. The hamiltonian of the model can be written in the form 

H = -Ji^criaj - J 2 ^2 0i °i ~ Jt^ViVjVkVi, (12) 

(ij) 

k I 

where Ji, J 2 and J 4 are the NN, NNN and plaquette couplings, respectively. This is a special 
case of the models studied in p^JT5|JT6|JT9|] . We shall first use the plaquette approximation 



of the CVM. Notice that this approximation has already been applied to the same model in 
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I|,[25|-p7|. In particular, Sanchez |^T| reported closed-form expressions for the equilibrium 



density matrices and the momentum space pair correlation function in the disordered phase. 



Moran-Lopez and coworkers p5| observed qualitatively the existence of a disorder locus in 



the phase diagram. Cirillo and coworkers [^7] calculated again the momentum space pair 



correlation function, and on this basis they determined the location of the disorder line. 



Their pair correlation function coincides with that by Sanchez except for a misprint p8 
and they obtained a disorder line which is only very close to the exact one, instead of 
coincident as it should be on the basis of the results of the present paper, because of an 
additional approximation. 

As a first step one can, at least at the numerical level, verify that the approximation is 
exact on the disorder variety using only the CVM. A simple way is to consider a hierarchy 



of approximations like the so-called C-series |29j , in which the maximal cluster is a rectangle 
made of 2 x L sites, with L > 2 (the plaquette approximation is the first element of the 
C-series). It is found, with extremely high precision, that a sequence of approximations in 
this series gives identical results on the disorder variety of the model. Inspection of the pair 
correlations shows that ([[J]) is also satisfied. 

On the other hand, using published results it requires only a long but straightforward 
calculation to check that on the (known |T2] , P^| , P^| ) disorder variety of the model one obtains 



the exact free energy. Looking at the pair correlations one also sees that the condition d = c 2 
(see (|TT1)) is satisfied on the variety of equation 

Uon exp(2J 4 )cosh(4J 2 ) +exp(-2J 2 ) 

COShl/Jl) = ; ; , (Id) 

V ' exp(2J 2 ) +exp(2J 4 ) ' V ; 



which is precisely the disorder variety of the model |l^J19|] . The free energy per site can be 
written as 

/ = - In [exp(- J 4 ) + exp( J 4 - 2J 2 )] , (14) 

and again coincides with the exact one, while the NN correlation is 

exp(-4J 2 ) - cosh(2Ji) 



sinh(2Ji) 

the NNN correlation d = c 2 and the plaquette correlation 

exp(4J 4 ) [1 - exp(8J 2 )] + 4exp(2J 2 ) [exp(2J 4 ) - exp(2J 2 ) 



(15) 



exp(4J 4 ) [1 - exp(8J 2 )] + 4exp(2J 2 ) [exp(2J 4 ) + exp(2J 2 )] ' ^ 



q = (<Ji<Jj<JkO-i) 



Finally, since all the pair correlations are given simply by (|TTf) we can easily calculate 
the momentum space correlation function, or structure factor. We first rewrite (pTTj) as 

(1^1 _|_ \y I \ 
J, where £ = — (lnc) -1 . After a Fourier transform one finds 

S(Px,P y ) = Si(p x )Si(p y ), where 

Qf\ sinh(l/Q 

cosh(l/£j — cosp 



It can be verified that the structure factors calculated by Sanchez |H] and (except for the 



misprint) Cirillo and coworkers [27] reduce to the above expression on the disorder line. 
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V. CONCLUSIONS 



I have shown that the CVM gives exact results on the disorder varieties of two- 
dimensional Ising-like models. In particular, I have considered the Ising model with NN, 
NNN and plaquette interactions on the square lattice, in the plaquette approximation of the 
CVM. In the disordered phase of the model, imposing the simple condition that the NNN 
pair correlation equals the square of the NN pair correlation, the CVM plaquette approxima- 
tion becomes exact and it is shown that this condition holds on the disorder variety, where 
the model can be solved in closed form. It is important to notice that, using the CVM, 
one can obtain any correlation function, due to the fact that the procedure of extending the 
local thermodynamic state is well-defined just on the disorder variety. Similar results can 
be obtained on the triangular lattice as well as other two-dimensional lattices. 
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FIGURES 

FIG. 1. The 3x3 square and the zig-zag chain 
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